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Abstract 

Background: A common assumption of microorganisms is that laboratory stocks will remain genetically and 
phenotypically constant over time, and across laboratories. It is becoming increasingly clear, however, that 
mutations can ruin strain integrity and drive the divergence or "domestication" of stocks. Since its discovery in 1960, 
a stock of Methylobacterium extorquens AM1 ("AM1") has remained in the lab, propagated across numerous growth 
and storage conditions, researchers, and facilities. To explore the extent to which this lineage has diverged, we 
compared our own "Modern" stock of AM1 to a sample archived at a culture stock center shortly after the strain's 
discovery. Stored as a lyophilized sample, we hypothesized that this Archival strain would better reflect the first-ever 
isolate of AM1 and reveal ways in which our Modern stock has changed through laboratory domestication or 
other means. 

Results: Using whole-genome re-seguencing, we identified some 29 mutations - including single nucleotide 
polymorphisms, small indels, the insertion of mobile elements, and the loss of roughly 36 kb of DNA - that arose 
in the laboratory-maintained Modern lineage. Contrary to our expectations, Modern was both slower and less fit 
than Archival across a variety of growth substrates, and showed no improvement during long-term growth and 
storage. Modern did, however, outperform Archival during growth on nutrient broth, and in resistance to rifamycin, 
which was selected for by researchers in the 1980s. Recapitulating selection for rifamycin resistance in replicate Archival 
populations showed that mutations to RNA polymerase B [rpoB] substantially decrease growth in the absence of 
antibiotic, offering an explanation for slower growth in Modem stocks. Given the large number of genomic 
changes arising from domestication (28), it is somewhat surprising that the single other mutation attributed to 
purposeful laboratory selection accounts for much of the phenotypic divergence between strains. 

Conclusions: These results highlight the surprising degree to which AM1 has diverged through a combination 
of unintended laboratory domestication and purposeful selection for rifamycin resistance. Instances of strain 
divergence are important, not only to ensure consistency of experimental results, but also to explore how 
microbes in the lab diverge from one another and from their wild counterparts. 
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Background 

To ensure that scientific results are both reproducible 
and consistent, a high level of integrity is required of ex- 
perimental methods and microbial stocks. One assump- 
tion is that stocks are constant over time, such that 
contemporary isolates of a given strain are genetically 
and physiologically identical across laboratories, and to 
stocks from many years ago. However, before the wide- 
spread use of deep freezers or lyophilization, storage of 
stocks using agar slants and other methods permitted 
growth and metabolism, albeit slowly, over long periods 
of time. And as long as stocks are metabolically active, 
mutations are likely to appear. If these mutations are 
beneficial, they can be enriched or fixed in stocks and 
subsequent sub-cultures via natural selection; or, alterna- 
tively, practices such as plate streaking and colony pick- 
ing could inadvertently propagate clones with neutral or 
even deleterious mutations due to random sampling or 
"genetic drift". Either way, such mutations are the bane 
of microbial stocks: they destroy the integrity of other- 
wise isogenic lines, and they become the raw material 
for processes such as selection and drift to facilitate evo- 
lutionary divergence in strains over time. 

The slow accrual of mutations often goes unnoticed in 
laboratory strains, but can result in considerable genomic 
and phenotypic differences both between independent la- 
boratory stocks, and between laboratory strains and their 
wild counterparts. This unintended mutational divergence 
is termed "laboratory domestication", and is particularly 
common in many (if not most) microorganisms isolated 
prior to the widespread use of modern storage methods 
(e.g., cryopreservation in deep freezers or lyophilization). 
For example, the domestication of Bacillus subtilis to the 
lab is associated with a loss of biofilms [1], swarming be- 
havior [2], and fruiting body formation (sporulation) [3]. 
Stocks of Salmonella enterica serovar Typhimurium [4,5] 
and Escherichia coli [6-8] archived for years or decades 
show a considerable amount of genetic diversity and novel 
phenotypes stemming from long-term growth and sur- 
vival during storage in agar stabs. Other controlled studies 
of microbes from the lab [9,10] or the wild [11-13] also 
show how growth, storage, and passaging procedures can 
readily lead to the diversification, divergence, and domes- 
tication of stocks. Outside of microbes, examples of do- 
mestication have been documented in stocks of the 
nematode worm, C. elegans [14,15], and in independent 
stocks of laboratory mice [16,17]. In all these examples, 
understanding the extent to which laboratory domestica- 
tion has occurred is important - not only for the 
standardization of experiments and the correct inter- 
pretation of results - but also because each instance of 
domestication is itself an interesting case-study of gen- 
omic and phenotypic divergence driven by a subtle and 
often cryptic set of evolutionary processes. 



Since the early 1960s, Methylobacterium extorquens 
AMI has emerged as the predominant model system for 
studies of bacterial one-carbon metabolism. As a faculta- 
tive methylotroph, M. extorquens AMI (herein referred 
to as "AMI") has the ability to grow using reduced one- 
carbon (Ci) compounds such as methanol and methyl- 
amine as the sole source of carbon and energy, as well 
as multi-carbon (multi-C) compounds such as succinate, 
pyruvate, and acetate [18,19]. The oxidation of methanol 
into biomass proceeds via the highly toxic intermediate, 
formaldehyde, and is complex, requiring over 100 en- 
zymes [20]. A sequenced genome [21], genetic tools 
[22-26], optimized growth conditions [27], metabolic 
models [28], and extensive knowledge of both Ci and 
multi-C metabolism [29] all make AMI the ideal organ- 
ism for studies of methylotrophy in the lab, as well as an 
emerging system for experimental evolution [30-33]. 
Aside from AMI, related methylobacteria are known for 
their roles in the plant microbiome [34-37], the biodeg- 
radation of toxic chemicals like chloromethane [38] and 
dichloromethane [39], and their potential for use in in- 
dustrial applications [40,41]. 

Unlike most organisms, the discovery and establishment 
of AMI as a model system was completely accidental. 
Around the year 1960, Dr. J.R. Quayle and colleagues at 
the University of Oxford were searching for a new organ- 
ism in which to study the oxidation and assimilation of Ci 
compounds, but discovered in their medium a "heavy, 
pink growth, presumably due to some airborne contamin- 
ant" [18,42]. After the contaminant was isolated, it was 
found to grow rapidly on a variety of both Ci and multi-C 
compounds. Shortly thereafter, a sample of AMI - then 
known as Pseudomonas AMI for "Airborne Methylo- 
troph #1" - was deposited to the National Collection of 
Industrial and Marine Bacteria (NCIMB, Aberdeen, 
Scotland), while a working stock remained in the lab. 
Over the years, this working stock was maintained and 
propagated between different researchers, laboratories, 
and growth and storage conditions, all the way up to our 
own laboratory's stock. These circumstances raise the 
question: to what extent has this AMI lineage diverged 
during its time in the lab? 

To address this question, we sought to compare to- 
day's AMI to its ancestor isolated circa 1960, or a close 
descendent of this ancestor. Closely related strains of 
M. extorquens differ significantly in their gene content 
and metabolic capabilities [43], making these a less 
than ideal comparison to determine the ancestral "wild" 
AMI state. Luckily, however, we realized that two major 
lineages of AMI were established circa 1960: the stock ar- 
chived circa 1961 at the culture stock center (herein re- 
ferred to as the "Archival" strain); and the working line of 
AMI that was propagated over many years from J.R. 
Quayle's lab, to Mary Lidstrom's group, to our own 
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laboratory stock of "Modern" AMI (Figure 1). We hy- 
pothesized that the Archival AMI - which underwent 
longer periods of lyophilized storage with fewer growth 
cycles in between - might better-reflect the ancestral state 
of AMI circa 1960, and offer an excellent reference with 
which to identify evolutionary divergence in the laboratory- 
maintained AMI lineage. We document here the surprising 
extent to which our Modern AMI has changed during its 
time in the lab using various assays of growth and fitness, 
long-term growth and storage, and whole-genome sequen- 
cing of the Archival AMI strain. We then offer a discussion 
of specific laboratory practices and evolutionary processes 
that may have led to such divergence. 

Results 

Whole-genome sequencing reveals extensive genomic 
divergence in the Modern AMI lineage 

To explore the extent to which AMI has diverged at the 
genomic level, we used whole-genome sequencing to 
compare the Archival genome to a previously sequenced 
Modern reference [21]. Illumina sequencing reads were 
analyzed both by mapping onto the Modern reference, 
and through de novo sequence assembly. For sites in 
which these strains differed, we compared the muta- 
tional state at each site (i.e., Archival or Modern) to 
other previously sequenced strains of M. extorquens to 
determine whether substitutions occurred in either the 
Modern or Archival lineage (Figure 1). While this ana- 
lysis cannot identify changes that occurred between the 
divergence of "wild" AMI from the Archival/Modern la- 
boratory ancestor, these mutations, if any, would only 
add to the ways in which AMI has evolved in the lab. 

Our results identified a sizeable number and variety of 
mutations that separate the Modern and Archival strains. 
We discovered 11 SNPs, 4 small indels (1-5 bp), the pro- 
liferation of 9 insertion sequence (mobile) elements, and 
some 36 kb of DNA found in 5 de novo assembled contigs 
that are present in the Archival strain but absent in 
Modern AMI (Table 1). For all but two of these muta- 
tions, excluding DNA loss, the Archival state is universally 



conserved with related M. extorquens strains to the exclu- 
sion of Modern AMI. Taken together, these results 
strongly suggest that this extensive genomic divergence oc- 
curred solely in the Modern AMI lineage, while the Arch- 
ival AMI has been largely preserved. 

Some Modern mutations are likely to have far-reaching 
physiological consequences. Certain mutations target 
highly pleiotropic genes such as rplj, which encodes 
the ribosomal subunit protein L10; rpoB, the beta sub- 
unit of RNA polymerase; and recG, the primary DNA 
helicase involved in recombination repair and other 
functions. The insertion of mobile elements might also 
have altered physiology in the evolution of Modern 
AMI, jumping into several putative protein-coding 
genes. Still other mutations resulted in the loss of a 
substantial amount of DNA from Archival to Modern 
AMI. Deleted regions were typically flanked by inser- 
tion sequences and other low complexity DNA, but in 
most instances we could infer the likely genomic loca- 
tion and content of these deletions (Table 1). Those 
genes deleted in Modern AMI are predicted to perform a 
variety of functions, and appear in some instances to be 
homologs of genes found only in distantly related mem- 
bers of the Methylobacterium genus. Further insights into 
the extensive loss of DNA in Modern AMI, as well as the 
functional and evolutionary consequences of Modern mu- 
tations, will require considerable future work. Still, these 
results clearly show that the Modern AMI lineage has 
undergone a substantial degree of evolutionary divergence 
and domestication over fifty years of growth and storage 
in the lab. 

Archival is faster and fitter under most standard 
growth conditions 

To begin to explore the phenotypic differences arising in 
Modern AMI, we compared the performance of Modern 
and Archival strains grown on standard growth sub- 
strates in "Hypho" minimal medium. Two primary met- 
rics were used to assess growth: specific growth rate, 
determined from the increase in optical density (OD 60 o) 
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Figure 1 Two distinct lineages of Methylobacterium extorquens AM1. Shortly after its discovery in 1960 (1), a sample of M. extorquens AM1 
("AMI") was deposited to a culture stock center for storage and distribution (2; Archival AM1). Many researchers, however, use instead a working 
stock of AM1 that was maintained over fifty years in the lab (3; Modern AM1), and was at one point selected for rifamycin resistance (Rif* 1 ) [45]. 
We hypothesized that these conditions may have fostered the accumulation of mutations and unintended evolutionary divergence in the 
Modern AM1 lineage, and sought to compare our Modern AM1 to the Archival strain. Dashes represent the accumulation of mutations in the 
Modern lineage. 
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Table 1 List of mutations derived in the Modem AMI lineage 





Single nucleotide polymorphisms 








Chromosome 


Position 


Mutation 


Annotation 


Gene/Locus 




Description 


META1 


Dl l,/UU 


c 
u - 


— > A 


intergenic 


META1_0578/META1_ 


0579 


Transcriptional accessory protein/ 
hypothetical protein 


META1 


2,050,899 


G- 




S359R (AGG -> AGO 


META1J984 




Putative catecholate siderophore 
receptor fill precursor (TonB- 
dependent receptor fiu) 


META1 


2,51 1,236 


C- 




L8004F (CIL — > I IL) 


Mb I A l_24lz 




Hypothetical protein 


MFTA1 


3,795,848 


c- 


->T 


L264F (CTG — > TTC) 


aldB 




Aldehyde dehydrogenase; 
chloroacetaldehyde dehydrogenase 


META1 


4,123,848 


c- 


-+T 


Intergenic (-1754/+2180) 


META l_4038/rffH 




Fragment of transposase related to 
IS701 family/glucose-1 -phosphate 
thymidylyltransferase 


MFTA1 


4,382,526 


c- 


-»G 


~F1 r r\~r l k r~ ^ A r~ ^\ 

T1 60T (ACC— > ACG) 


A ACTA T A ~)m 

Mt I A l_4292 




Plasmid replication protein RepA 


METTAI 


4,494,119 


A- 


^G 


N55D (/IAC — > GAC) 


rpU 




50S ribosomal subunit protein L10 


MFTA1 


4,496,733 


A- 


^G 


Q521R(CAG->CGG) 


rpoB 




RNA polymerase, beta subunit 


META1 


4,498,413 


A- 


->G 


Q1081R(G4G^CGG) 


rpoB 




RNA polymerase, beta subunit 


META1 


4,665,003 


C- 


-> A 


L300M (CTG->ATG) 


recG 




DNA helicase, ATP-dependent 
resolution of Holliday junctions, 
branch migration 


META1 


5,187,175 


G- 


-> A 


A260T(GCC->-A;Q 


METAl_5044 




Putative o-succinylbenzoate- 
CoA igase 




Small insertions & deletions 










Chromosome 


Position 


Mutation 


Annotation 


Gene/Locus 




Description 


META1 


1,083,921 


+1 bp 


Intergenic (-299/+188) 


METAl_l04l/METAl_ 


1 042 


Putative CoxB/conserved 
hypothetica protein 


META1 


2,012,135 


+5 bp 


Ceding (126-130/228 nt) 


META1J939 




Hypothetical protein 


META1 


5,018,430 


A1 bp 


Intergenic (-120/+218) 


META 1_4900/META l_490l 


Putative hydrolase of beta-lactamase 
superfamily/conserved hypothetical 
protein DUF949 


META2 


821,910 


+1 


bp 


Intergenic (+133/-26) 


META2_0S63/META2_ 


0864 


AAA superfamily ATPase/ 
hypothetica protein 




IS elements 














Chromosome 


Position 


Mutation 


Annotation 


Gene/Locus 




Description 


META1 


772,350 


+1408 bp 


Gain !SMex3; coding 


METAL0742-METAI 


0743 


META 1_0742, META1_0743 


META1 


929,023 


+1390 bp 


Gain SMex14; intergenic 


META1_0895-META1_ 


0898 


Transposase of ISMexl, IS3 family 
(ORF 1)/transposase of ISMexl 4, 
IS256 family 


META1 


3,730,805 


+1 


408 bp 


Gain SMex3; coding 


META1_3592-METAI 


_CDS3732187D 


META 1_3592, META 1_3593, 
META1_CDS3732187D 


META1 


4,143,329 


+1620 bp 


Gain ISMex4; intergenic 


META1_4059-META1_ 


4061 


Hypothetical protein/transposase of 
ISMex4, IS1 380 family 


i\ A CT A 1 

IVlt 1 A 1 


4,149,803 


+1408 bp 


Gain ISMex3; coding 


META 1_4069 




IVlt IA l_4Uuy 


META1 


4,702,223 


+1620 bp 


Gain ISMex4; intergenic 


META1_4586-META1_ 


CDS4704205D 


Transposase of ISMex4, 151380 
family/hypothetica protein; 
RMQ08497 


META1 


4,909,262 


+1620 bp 


Gain ISMex4; coding 


META 1_480 1-META1_ 


.4803 


Hypothetical protein///transposase 
of ISMex4, IS1 380 family/conserved 
hypothetical protein 


META2 


426,304 


+1205 bp 


Gain ISMexl; coding 


META2_0472-MEJA2_0475 


Transposase of ISMexl, IS3 family 
(ORF 1)/conserved hypothetical 
protein 


META2 


1,153,967 


+1620 bp 


Gain ISMex4; intergenic 


META2J243-MEJA2_ 


1245 


Transposase of ISMexl 0, ISL3 family/ 
transposase of ISMex4, IS1380 family 
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Table 1 List of mutations derived in the Modem AMI lineage (Continued) 



Unmapped Archival AMI DNA 


Size (bp) 


Location(s) 


Composition 


11908 


? 


Porin protein, transcriptional regulator (AraC) protein, conjugative relaxase domain protein, sodium/hydrogen 
exchanger, TraG homolog 


8419 


META1 _4345/META2_01 37 


TonB-dependent receptor/siderophore receptor protein, hypothetical proteins 


8000 


META1J083 


Sodium/calcium exchanger, hypothetical proteins 


5207 


META2_0137 


Cold shock protein A {cspA), metallophosphoesterase, plasmid stabilization system, addiction module antidote 
protein, cobyrinic acid ac-diamide synthase, stability/partitioning determinant, hypothetical proteins 


2423 


p2META_0017 


Oxidoreductase molybdopterin binding protein, sulfiterytochrome c oxidase subunit B, hypothetical proteins 


35957 


Total bp 





Mutations were identified by comparing the Archival genome to a previously sequenced Modern reference [21]. By comparing the mutational state {Archival or Modern) at 
each site to other previously sequenced strains of M. extorquens (see Methods), all but two mutations can be unambiguously traced as having occurred on the branch from the 
Ancestral to Modern AMI (Figure 1). The effect of nonsynonymous mutations on coding regions are highlighted in bold and italic. 



for a particular strain/condition using a custom-built, 
high-throughput microbial culturing system and analysis 
software [27,44]; and relative fitness, which analyzes 
performance over multiple phases of bacterial growth 
(i.e. lag, exponential, and stationary) using a head-to- 
head competition of strains in co-culture [30]. 

Our initial hypothesis was that Modern AMI would out- 
perform the Archival strain, owing to an increased likeli- 
hood of mutations in Modern AMI that could facilitate 
adaptation to laboratory conditions. However, contrary to 
our expectations, we found that the Archival strain was 
both faster and fitter under most conditions tested. The 
Archival strain was considerably faster growing on the Ci 
compounds methanol (42%) and methylamine (12%), as 



well as the multi-C substrate succinate (52%; Figure 2A 
and B). In head-to-head fitness competitions, the Archival 
strain showed roughly a 30% advantage across the sub- 
strates tested (Figure 2C), suggesting that analyses of 
growth rate and fitness are not entirely correlated. None- 
theless, these results show that the Archival strain outper- 
forms Modern under most standard growth conditions, 
and suggest that the Modern lineage became slower and 
less fit during its time in the lab. 

Modern outperforms the Archival AMI on nutrient broth 

We next sought to compare Modern and Archival AMI 
using less traditional growth substrates. In contrast to 
"Hypho" minimal medium, nutrient broth (NB) is a rich 




20 

Time (hrs) 



M M Ma S 

Figure 2 Archival outperforms Modern AM1 under standard growth conditions. A) Representative growth curves showing the increase in ODan 
over time of Modern (black circles) and Archival AM1 (gray squares) cultured using 3.5 mM succinate in 48-well plates. B) Growth rates calculated from 
the exponential phase of cultures grown on methanol (M), methylamine (Ma), or succinate (S) as a carbon source. Significant growth differences between 
Modern and Archival were calculated using a two-tailed, unpaired t test, and are marked by single (p < 0.05) and double asterisks (p < 0.01). C) Fitness of 
Archival AM1 measured via a head-to-head competition mixed in co-culture with a fluorescently labeled Modern reference. A control 
growth consisted of unlabeled Modern (black) versus the fluorescent Modern reference grown on M. All other bars (gray) show Archival 
fitness relative to Modern grown M, Ma, and S. Values are the mean plus SEM of growth rates or fitness values calculated from three or more 
biological replicates (see Methods). 
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medium composed of partially digested proteins (pep- 
tone), but its exact nutritional components are largely 
undefined. Today, NB is rarely used as a growth sub- 
strate for Methylobacterium, and is used almost exclu- 
sively to support growth in co-cultures with E. coli 
during conjugal matings for genetic manipulation. It is 
possible, however, that Modern AMI was more fre- 
quently cultured on NB in the past and may have 
adapted to growth on this medium. To test this hypoth- 
esis, we assessed the growth and fitness of Modern and 
Archival AMI cultured on NB. 

We found over multiple different growth experiments 
that Modern AMI consistently outperformed Archival 
on NB, although both strains grew more poorly than on 
Hypho medium. Modern displayed a clear growth ad- 
vantage on NB in 48-well plates (Figure 3A), which 
could be traced in part to its relative insensitivity to per- 
turbations arising during the later stages of NB growth. 
While monitoring the increase in OD 60 o over time, we 
observed that NB cultures tended to slow down during 
the later stages of exponential growth, and that the 
Archival strain was hindered to a greater extent than 
Modern. This suggested that Modern may have adapted 
to yet unknown components of growth in NB; however, 
due to this decrease during late exponential phase, we 
were unable to accurately assess differences in specific 
growth rate, and sought instead to quantify performance 
using a head-to-head competition of Modern and Arch- 
ival co-cultures. 

The performance of Archival AMI co-cultured with a 
fluorescently labeled Modern reference was monitored 
over the course of several days during growth in NB 
flasks. Co-cultures were sampled periodically to monitor 
changes in the ratio of nonfluorescent (Archival) to 
fluorescent (Modern) cells using flow cytometry. These re- 
sults suggest that the Archival AMI holds an early growth 
advantage in NB that quickly decreases from 1 to 4 days 
until the strains reach stationary phase (Figure 3B). As a 
control, a co-culture of nonfluorescent Modern mixed with 
the same fluorescent Modern reference showed little 
change over the course of the experiment. Using the ratio 
of nonfluorescent to fluorescent cells at the start and the 
end of one growth cycle, we can calculate the fitness of the 
Archival AMI relative to Modern assuming a 64-fold (2 ) 
increase in the population (see Methods), and find that 
Archival is only 61% as fit as Modern during growth on 
NB. Overall, the improved performance of Modern on NB 
is suggestive of adaption, either specifically to this medium, 
or to yet unknown but similar growth conditions. 

Archival and Modern are similar in terms of long-term 
growth and storage 

Another dimension in which AMI may have adapted to 
life in the lab is through improved performance during 




0 I i I i I i I i | 

0 12 3 4 

Days in culture 

Figure 3 Modern outperforms Archival AMI when grown on 

nutrient broth. A) Representative growth curves of Modern (black 

circles) and Archival (gray squares) AM1 grown on nutrient broth 

(NB). Note that growth - particularly for the Archival strain - slows 

considerably during late exponential phase, signifying density-dependent 

growth inhibition. B) Change in the proportion of either Modern or 

Archival AM1 mixed in co-culture with a fluorescently labeled 

Modern reference as measured by flow cytometry. Values represent 

the mean plus SEM of at least three biological replicates grown in 

48-well plates (A) or flasks (B). 
\ ) 

long-term growth and storage. We compared the Mod- 
ern and Archival strains grown for extended periods 
both in flasks and on plates by creating co-cultures of 
each strain with a fluorescently labeled Modern refer- 
ence, and monitoring the change in fluorescent to non- 
fluorescent cells over time using flow cytometry. Growth 
in flasks was performed over 14 days with continual 
shaking using succinate as a growth substrate, while 
growth on methylamine plates was carried out at 30°C 
for 4 days, followed by up to 60 days at 4°C to simulate 
long-term storage. 

Here, Modern and Archival AMI were roughly equiva- 
lent in terms of growth and survival during stationary 
phase. In flasks, the ratio of unlabeled Archival cells 
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to labeled Modern cells was steady for 12 days after 
the initial 2 days of exponential growth (Figure 4A), 
while the unlabeled Modern control remained un- 
changed for the duration of the experiment. Archival 
held a similar advantage during storage at 4°C on agar 
plates for up to 60 days (Figure 4B). Collectively these 
results suggest that, at least under these conditions 
tested, the extent to which laboratory domestication 
improved the long-term growth and survival of AMI 
is limited, and that the major difference between these 
strains lies in decreased exponential phase growth in 
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Figure 4 Equivalence of AMI strains during tests of long-term 
growth and survival. Co-cultures were created by mixing either 
Modern (black circles) or Archival AM1 (gray squares) with a fluorescently 
labeled Modern reference, and the change in unlabeled versus fluorescent 
cells was monitored over time using flow cytometry. A) In continually 
shaken flasks with succinate, the Archival strain increased in frequency 
over the first two days of growth and maintained this advantage over 
Modern over time. B) Similarly, Archival increased in frequency during four 
days of growth on methylamine agar plates (not shown), and maintained 
this frequency during long-term storage at 4°C. Values represent the 
mean plus SEM of the percent unlabeled cells measured in three 
replicate co-cultures. 



Modern under standard conditions, and slightly in- 
creased growth on NB. 

Growth of AM1 is significantly hampered by selection for 
rifamycin resistance 

To explain the decreased performance of Modern AMI 
under most growth conditions, we returned to the gen- 
omic changes identified in this lineage. All but one of 
the 29 mutations arising in Modern AMI can be attrib- 
uted to unintended laboratory domestication; this single 
exception, however, was central to the development of 
Modern AMI. In 1984, Fulton and colleagues selected 
for rifamycin resistance (Rif ) to facilitate genetic manip- 
ulations in AMI using conjugative, tri-parental matings 
[45]. Across numerous systems, mutations conferring 
Rif most often occur in the beta subunit of RNA poly- 
merase - encoded by rpoB - and typically give rise to 
fitness tradeoffs between survival in the presence of anti- 
biotics versus decreased growth in their absence [46]. 
Compared to the Archival strain, Modern AMI has two 
mutations to the rpoB locus (Table 1), one of which 
(Q521R) falls in a region that confers Rif 8 in a variety of 
other organisms [47], while the effect of the other 
(Q1081R) is yet unknown. At both positions the amino 
acid state of the Archival strain is universally shared 
with other non-AMl strains of M. extorquens, suggest- 
ing that these mutations arose exclusively in the Modern 
lineage. Thus, these mutations to rpoB, particularly 
Q521R, are strong candidates for decreased overall 
performance in Modern AMI: offering resistance in 
the presence of Rif but slower growth in its absence. 

To explore the potential costs associated with anti- 
biotic resistance, we recapitulated the evolution of Rif R 
in independent replicate cultures of the Archival strain. 
Out of 36 independent populations grown to saturation, 
only 7 produced a small number of spontaneous, resist- 
ant colonies when plated on Rif agar plates. Each inde- 
pendent population was streaked to purity, analyzed in 
terms of growth rate on succinate (with no antibiotic), 
and sequenced along with the Modern and Archival 
strains at the rpoB locus. 

Independent experiments recapitulating Rif R in the 
Archival strain all selected for mutations to rpoB, and all 
resulted in decreased growth in the absence of antibiotic. 
Upon sequencing rpoB from the 7 newly evolved Arch- 
ival isolates, we found that Rif R was always associated 
with mutations to rpoB. One of these strains (CM4022) 
acquired the exact same nonsynonymous change that 
occurred in Modern evolution, corroborating our hy- 
pothesis that the Q521R mutation (from Archival to 
Modern) is causal in conferring Rif R to Modern AMI. 
Indeed, all mutations observed fall within a region of 
rpoB that is commonly mutated to confer Rif R across a 
variety bacteria [47]. 
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The spectrum of rpoB mutations across Rif K Archival 
isolates displayed highly variable effects on growth rate 
in the absence of antibiotic. Compared to their Archival 
ancestor, several Rif isolates show very little decrease in 
performance when grown on succinate, while other iso- 
lates slow to near Modern levels, and still others grow 
substantially worse than Modern (Figure 5B). Interest- 
ingly, the CM4022 isolate that perfectly recapitulated the 
change from Archival to Modern AMI (Q521R) was 
slightly faster than Modern, suggesting perhaps that 
other mutations further hamper growth in the Modern 
lineage. We note, however, that a direct comparison of 
this strain is difficult given that many other mutations 
and mutational interactions were likely present in Mod- 
ern AMI during the original selection for Rif R . Never- 
theless, these results demonstrate that selection for Rif R 
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Figure 5 Mutations associated with rifamycin resistance hinder 
AM1 growth. A) Spectrum of mutations to the RNA polymerase 
beta subunit (RpoB) during past and current selection for rifamycin 
resistance (Rif 3 ). Modern AMI was selected for Rif" in 1984 [45] and 
acquired two mutations to RpoB: Q521R, and Q1081R (denoted by 
an asterisk). By recapitulating selection for Rif" in replicate Archival 
populations, we identified a number of other RpoB mutations 
putatively associated with Rif^. B) The effect of Rif mutations on 
growth rate in the absence of antibiotic. Values represent the mean 
plus SEM of four biological replicates grown in 48-well plates with 
succinate. Strains that were significantly slower than Archival are 
marked with asterisks (p < 0.01, one-way ANOVA with Tukey 
post-hoc test). 



can substantially hinder growth of the Archival strain 
in the absence of antibiotic, and that this single 
researcher-imposed event - not laboratory domestica- 
tion - is the major mechanism by which Modern AMI 
became slower growing in the lab. 

Discussion 

This work highlights the surprising extent to which 
M. extorquens AMI has inadvertently diverged during 
fifty years of growth and storage in the lab. Compared 
to an Archival AMI strain, our Modern AMI stock 
was slower growing under most standard laboratory 
conditions, and one mechanism to account for this 
decrease is in past selection for rifamycin resistance. 
Indeed, the recapitulation of Rif R in independently 
evolved populations of the Archival strain frequently, 
but not always, led to a tradeoff between survival in 
the presence of antibiotic and decreased growth in its 
absence. Upon sequencing the Archival genome, we iden- 
tified some 29 mutations that have accrued in the Modern 
AMI lineage, including a number of single nucleotide 
polymorphisms, small insertions and deletions, the prolif- 
eration of mobile elements, and the loss of some 36 kb 
of DNA. Though the full impact of these mutations for 
improving growth on NB or other conditions has yet to be 
revealed, it is clear that Modern AMI has diverged sub- 
stantially through laboratory domestication and changes 
stemming from selection for antibiotic resistance. 

At first glance, it was not immediately clear why Mod- 
ern AMI had become slower and less fit under standard 
growth conditions. Our initial hypothesis was that Mod- 
ern would outperform the Archival strain due to the ac- 
quisition of mutations that optimized growth or survival 
in the laboratory, but our results suggested the contrary. 
Alternative hypotheses to account for decreased per- 
formance in Modern might include the chance acquisi- 
tion of one or more mutations that are deleterious for 
growth, or that selection for mutations that improve 
growth under certain laboratory conditions generated 
tradeoffs in others. In the latter case, a strong candidate 
emerged in tradeoffs generated by selection for antibiotic 
resistance. Antibiotic resistance in organisms often car- 
ries with it a fitness cost in the absence of antibiotic, and 
this has been documented in related strains to our Mod- 
ern AMI [48]. We discovered two mutations to rpoB in 
our Modern strain: one that is hypothesized to confer 
resistance (Q521R), while the effect of the other is yet 
unknown but might serve a compensatory function [49]. 
Reselecting for Rif R in 7 independent populations of 
Archival AMI resulted in a variety of mutations to rpoB 
with variable effects on growth in the absence of anti- 
biotic. One strain in particular (CM4022) perfectly reca- 
pitulated the Q521R mutation in Modern AMI, and 
yet showed slightly faster growth than Modern despite 
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lacking any other known changes. While a direct com- 
parison of this strain and the original selection for Rif R 
in Modern AMI is complicated by the presence of other 
mutations in the genetic background of the latter, these 
results strongly suggest that selection for Rif R was the 
major factor affecting growth in the Modern strain. 
Given the large number of genomic changes identified 
in the Modern lineage (n = 29), it is somewhat surprising 
that the only mutation attributed to purposeful labora- 
tory selection, and not domestication, accounts for much 
of the phenotypic divergence observed between the 
Archival and Modern strains. 

Though the exact circumstances under which Modern 
AMI evolved cannot in most cases be ascertained, we 
can at least hypothesize as to the general factors that 
might have played a role. Prior to the widespread use 
of -80°C freezers, long-term storage of AMI cultures was 
accomplished using Hypho agar slants with methanol or 
methylamine [18,50], most often under refrigeration. A 
study of Salmonella [4,5] and E. coli [6-8] archived for 
decades under similar conditions revealed a number of 
genomic and physiological changes that aid in survival 
during long-term stationary phase. These conditions 
often select for mutants with a "growth advantage in 
stationary phase" (GASP) phenotype associated with 
increased catabolism of amino acids and other small 
organic compounds, as well as the ability to outcompete 
"younger" cultures [51]. Although a GASP-like phenotype 
was not observed in the two test environments with which 
we compared Modern versus Archival AMI, further 
experiments might reveal conditions in which Modern 
AMI has adapted to other facets of laboratory growth 
or storage. One particularly interesting direction is in 
the advantage of Modern on NB, which may have 
evolved specifically to this medium or merely correlate 
with improvements in other conditions, such as the 
ability to survive long-term growth with limiting nutrients. 
Along these lines, we analyzed the growth of Rif* strains 
(CM4020-26) on NB and found that they performed 
equally as well, if not worse, than their Archival parent 
strain (data not shown). This suggests that improved 
growth of Modern on NB is due to one or more mutations 
that arose during laboratory domestication, not through 
the rpoB mutation associated with Rif*. 

Between long periods of storage, Modern AMI may 
have also adapted to yet unknown growth conditions. 
The specific components of growth media can act as a 
strong selective pressure in microbial cultures [52-54], 
and over many years Modern AMI experienced both 
numerous variants of minimal media, as well as occa- 
sional growths in rich media like NB. Large populations 
of microorganisms competing for limiting resources can 
create strong selective pressures for increased growth 
rate, which may have accrued in Modern AMI only to 



be partly nullified due to later tradeoffs with Rif R . It is 
also possible that some mutations were fixed not 
through selection for improved performance in labora- 
tory conditions, but rather via genetic drift stemming 
from practices that result in extreme population bot- 
tlenecks, such as colony picking. Given the complex 
and overall vague history of Modern AMI in the lab, 
reconstructing the exact mechanisms that drove its 
evolution might be best accomplished by studying 
other "offshoots" of the laboratory-maintained Modern 
lineage, and by characterizing the selective effect of 
Modern mutations in the Archival genetic background 
using allelic exchange. 

Conclusions 

The laboratory environment affords researchers with a 
great degree of control over experimental variables: from 
reagents and protocols, to the genotype of their model 
microbe, and the environment in which it lives in. Ex- 
amples of laboratory domestication, however, highlight 
the difficulty of maintaining high quality microbial stocks. 
Mutations can jeopardize the integrity of microbial stocks 
and, given time, lead to spurious and inconsistent results 
stemming from the evolutionary divergence of strains. 
Even purposeful laboratory experiments intended merely 
to select for antibiotic resistance or otherwise alter the 
genetic background of strains can have unanticipated 
pleiotropic effects. Thus, extra care should be taken to en- 
sure that experimental findings from strains are reprodu- 
cible and consistent over time. For those stocks in which 
divergence has already occurred, these microbes offer a 
unique opportunity to explore genetic and phenotypic 
changes resulting from complex evolutionary processes at 
work in the lab. 

Methods 

Bacterial strains & growth conditions 

Strains relevant to this study were as follows. Our Mod- 
ern AMI strain was derived from a pink, "wildtype" M. 
extorquens AMI (CM501) described elsewhere [23]. A 
sample of lyophilized Archival AMI (renamed CM3944) 
was acquired from the National Center of Industrial 
Food and Marine Bacteria (NCIMB #9133, Aberdeen, 
Scotland), grown to saturation, and frozen at -80°C in 
8% DMSO. To limit cell clumping and reduce noise in 
OD 60 o measurements during analyses of growth rate, we 
utilized a strain of Modern AMI, CM2720, that lacks 
genes for cellulose biosynthesis [27] without a loss of 
growth or fitness. For competitions in co-culture, a fluo- 
rescently labeled Modern reference (CM1175) was con- 
structed by placing the red fluorescent protein mCherry 
under control of a constitutive P tac promoter at the 
katA locus [30]. Isolates from each of seven Archival 
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populations selected for rifamycin resistance (described 
below) were numbered CM4020 through CM4026. 

Standard growth conditions utilized a modified version 
of Hypho minimal medium consisting of: 100 mL phos- 
phate salts solution (25.3 g of K 2 HP0 4 plus 22.5 g 
Na 2 HP0 4 in 1 L deionized water), 100 mL sulfate salts 
solution (5 g of (NH 4 ) 2 S0 4 and 2 g of MgS0 4 • 7 H 2 0 
in 1 L deionized water), 799 mL of deionized water, and 
1 mL of trace metal solution [55]. All components were 
autoclaved separately before mixing under sterile condi- 
tions. Carbon sources added just prior to inoculation in 
liquid minimal media consisted of 20 mM methanol, 
3.5 mM sodium succinate, or 20 mM methylamine 
hydrochloride. Growths in 48-well microtiter plates con- 
sisted of Hypho medium plus the appropriate carbon 
source to a volume of 640 uL. Agar plates consisted of 
growth medium plus either 125 mM succinate or 
100 mM methylamine and were autoclaved with 1.6% w/v 
agar. Difco nutrient broth (Becton, Dickson, and Com- 
pany, Franklin Lakes, NJ) was prepared according to the 
manufacturer's guidelines. 

All growth regimes consisted of three phases consist- 
ing of inoculation, acclimation, and experimentation 
growths. All strains were stored in vials at -80°C in 8% 
DMSO; growths were initiated by transferring 10 uL 
freezer stock into 10 mL of Hypho medium with metha- 
nol. Upon reaching stationary phase (~2 days at 30°C 
with shaking), cultures were transferred 1:64 into the ap- 
propriate medium and vessel to be tested, allowed to 
reach saturation in this acclimation phase, and diluted 
1:64 again into fresh medium for the measured (experi- 
mental) growth. 

Measurements of specific growth rate and relative fitness 

The increase in OD 60 o for strains grown in 48-well mi- 
crotiter plates was measured using an automated, ro- 
botic culturing and monitoring system [27,44]. The 
specific growth rate of cultures was calculated from the 
log-linear growth phase using custom-designed growth 
analysis software [27]. Growth rates reported for each 
strain and condition are the mean plus SEM calculated 
from triplicate biological replicates, unless otherwise noted. 
Exogenous cellulase enzyme from Aspergillus niger (Sigma- 
Aldrich, St. Louis, MO) was added to the medium to a final 
concentration of 0.1 mg/mL to further minimize cell 
clumping and facilitate more accurate measurements of 
OD 60 o over time (SMC, unpublished). 

Fitness measurements - which encapsulate the lag, ex- 
ponential, and stationary phases of growth - were per- 
formed using a head-to-head competition of strains 
grown in co-culture [30]. Modern and Archival AMI 
were competed against a fluorescently labeled Modern 
strain (CM1175) expressing mCherry [30]. Co-cultures 
were prepared by mixing test strains with the fluorescent 



Modern reference in roughly equal optical densities. A 
sample of this co-culture prior to competition was di- 
luted in 8% DMSO and stored at -80°C; the rest was di- 
luted 1:64 into 640 uL Hypho medium plus carbon in 
48-well microtiter plates and incubated with shaking at 
30°C for 1 growth cycle. A sample of co-culture after 
competition was frozen for later analysis using flow 
cytometry. 

The ratio of labeled to unlabeled cells in co-cultures 
before (R0) and after (Rl) competition was measured 
using a BD LSR Fortessa flow cytometer with an HTS at- 
tachment (BD Biosciences, San Jose, CA). Both forward 
and side scatter settings were set to 300 V to account 
for the small size of bacterial cells [56], and the flow- 
rate was adjusted to the lowest setting to more accur- 
ately identify events (cells) in dilute co-cultures. The 
fitness (W) of strains relative to the labeled Modern 
reference was calculated using the following formula, 
which assumes a 64-fold expansion of cells following 6 
doublings per growth cycle: 



Analysis of long-term growth and storage 

To assess the ability of strains to withstand extended pe- 
riods of stationary phase growth in flasks, Modern and 
Archival AMI were mixed in co-culture with the fluorescent 
Modern reference, and the ratio of labeled to unlabeled cells 
was measured periodically using flow cytometry. Flasks pos- 
sessing succinate were sealed to limit evaporation, and at 
the conclusion of the experiment co-cultures were streaked 
onto nutrient agar plates to check for contamination. 
Growth and survival on plates was measured using a similar 
experimental design: co-cultures were spread onto Hypho 
plus methylamine plates, grown for 4 days at 30°C, and then 
stored at 4°C for up to 60 days with periodic sampling. 

Recapitulation of Rif resistance in AMI 

The re-evolution of Rif resistance in Archival AMI was 
performed using 36 replicate lineages grown from single 
colonies in 48-well plates. After 2 days growth in liquid 
Hypho medium plus succinate, cultures were plated 
without dilution onto Hypho agar plus succinate plates 
containing 50 ug/mL Rif. After 5 days of growth, Rif R 
colonies were obtained from 7/36 cultures and streaked 
twice more to ensure the purity of clonal isolates 
(CM4020-4026). PCR amplification plus sequencing was 
used to assess mutations to the rpoB locus and the 
growth rate of strains on succinate was determined as 
described above. 
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Whole-genome sequencing of Archival AMI 

Preparation of Archival genomic DNA was performed as 
described in [57]. Briefly, an Archival cell pellet was 
lysed using bead beating (MP Biomedicals, Solon, OH) 
and digested using heat, SDS and Triton-XlOO deter- 
gents, proteinase K, and RNase A. Cell debris was pre- 
cipitated with CTAB/NaCl, removed using a phenol/ 
chloroform/isoamyl alcohol extraction, and genomic 
DNA recovered using an isopropanol precipitation. Illu- 
mina HiSeq2000 sequencing was performed out-of- 
house using GENEWIZ, Inc (South Plainfield, NJ). A 
total of 23,660,656 reads were generated, and 98.2% of 
these were mapped against the Modern AMI reference 
[21] using breseq v0.21 [58] with Bowtie 2 version 2.0.0- 
beta7 [59]. To identify DNA lost in Modern AMI evolu- 
tion, reads that were unmapped using breseq were 
analyzed using ABySS vl.3.4-kmer-96 [60] for de novo 
assembly of contigs,. A full de novo assembly of the 
Archival genome was also performed in ABySS for com- 
parison against breseq. Contigs from ABySS were further 
assembled and curated using contig assembly and 
BLAST in Geneious [61]. To determine whether muta- 
tional differences occurred in either the Modern or 
Archival lineage, we compared the mutational state 
(Archival or Modern) at each site to other previously 
sequenced strains of M. extorquens using both BLAST 
and analyses of conserved gene synteny using the 
Methyloscope project in MicroScope [62]. Mutations were 
identified as "derived" in the Modern lineage if the Archival 
state was consistency shared across related strains, to the ex- 
clusion of Modern AMI. Related strains included M. extor- 
quens PA1 (GenBank Assembly ID: GCA_000018845.1), 
CM4 (GCA_000021845.1), DM4 (GCA_000083545.1), and 
DSM 13060 (GCA_000243435.2). For all loci but two, the 
Archival state was universally conserved; the only exceptions 
were for genes with little to no homology outside of Modern 
and Archival AMI. 

Availability of supporting data 

Whole-genome sequencing data of the Archival AMI 
strain has been deposited to the NCBI Sequencing Read 
Archive (SRR 1046370). 
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